Spin ice on the trillium lattice studied by Monte Carlo calculations 
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We study a local ferromagnetic Ising model for classical spins on the trillium lattice. The ground 
state of this model features two spins out(/in) and one spin in(/out) on each triangle, and leads to 
a macroscopic ground state degeneracy. Our Monte Carlo simulations find a ground state entropy 
intermediate to that of spin ice on the kagome and pyrochlore lattices, suggesting that trillium spin 
ice is highly frustrated. To motivate the search for trillium spin ice, we calculate the magnetic 
susceptibility and structure factor. We note the qualitative resemblance of the susceptibility to 
previously published work on EuPtSi, which features local moments on the trillium lattice. 

PACS numbers: 75.10.Hk, 75.40.Cx, 75.40.Mg 



I. INTRODUCTION 

Recently there has been great interest in the possibil- 
ity that spin ice materials may provide a route to the 
realization of deconfined magnetic monopoles in a fully 
three dimensional (3D) correlated spin system. While 
such proposalsii^ and their experimental support^i^ have 
focused on materials based on the pyrochlore structure, 
a growing stable of lattice structures, and accompany- 
ing spin ice materials, is developing which may provide 
complementary routes to the same fundamental physics. 
In particular the existence of (two-dimensional) kagome 
spin ice was predicted in 2002 by Wills et al.—, extended 
to dipolar spin ice by Chcrn et al^, and recently real- 
ized with dipolar mcsoscopic interactions between litho- 
graphically etched ferromagnetic islands^. The existence 
of (3D) hyperkagome dipolar spin ice has recently been 
proposed^ as a pure limit of "diluted" spin ice^. 

Spin ice on the pyrochlore lattice is one of the most well 
understood examples of a geometrically frustrated (GF) 
system both theoreticalljiiSrJ^ and experimentallyi^"— . 
GF materials have particularly interesting properties 
which arise because the symmetries of the triangle- or 
tetrahedral-based lattice structures lead to interactions 
between charges and/or spins which cannot simultane- 
ously be uniquely minimized. Gcnerically one finds many 
equal energy ground state configurations among which 
selection of a long range ordered ground state manifold 
proceeds via higher order processes such as order by dis- 
order produced by finite temperature or quantum fluctu- 
ations. One of the earliest recognized frustrated materi- 
als was water, which was found to retain a finite entropy 
to very low temperatures by Giauque and Stout^ and ex- 
plained in terms of proton disorder by PaulingSS. A short 
14 years later it was recognized that macroscopic ground 
state entropies could arise in models of magnetism on GF 
lattices when Wannier— showed that antiferromagnetic 
(AFM) Ising spins on the triangular lattice had a ground 
state entropy per spin of 0.338314. By the mid 1980s 
it was recognized that a number of Ising AFMs exhibit 
finite ground state entropies^^ and the list of GF lattices 
exhibiting finite ground state entropies has continued to 
grow2^, now including some quantum spin models in ad- 



dition to the spin ice materials motivating this work. 

The term "spin-ice" was coined by Harris et alr^ to de- 
scribe the physics resulting from a ferromagnetic (FM) 
coupling between nearest neighboring spins on the py- 
rochlore lattice when these local moments experienced 
an Ising anisotropy directed along the local < 111 > di- 
rections (which point toward the centre of each tetrahe- 
dron). Following earlier work by Anderson^!, Harris et 
al. recognized that the ground state spin directions of 
such a local FM Ising model had a one-to-one correspon- 
dence with the displacements of hydrogen atoms from 
a corner-shared tetrahedral lattice in cubic water iceS^. 
This correspondence enabled Harris et al. to find a sim- 
ple approximation to the residual entropy of Ho2Ti207 
based on Pauling's estimate^S for the residual entropy of 
water ice made almost 60 years earlier. 

On the pyrochlore lattice the Pauling estimate for the 
ground state entropy per spin of spin ice is ^ = ^ ln(|). 
Pauling's estimate results from considering the probabil- 
ities of different ground state occurrences to be uncor- 
rclated. To obtain this result, one notes that of the 2^ 
spin possibilities for N Ising spins, only six of the six- 
teen possible configurations on each tetrahedron belong 
to the ground state with two spins in and two spins out. 
This reduces the degeneracy of the ground state by a 
factor of {jq)^^ where A^a = y number of tetra- 

hedra composing the lattice, such that for uncorrelated 
spins one expects a ground state degeneracy of 2''^(|)^. 
Since Harris et aZ.'s early work it has been realized that 
dipolar interactions also play a strong role in the experi- 
mental realization of spin-ice materials on the pyrochlore 
latticeiSiii, with the possibility of creating an effectively 
FM coupling between nearest neighbor spins even when 
the Ising coupling has an antiferromagnetic sign. 

Wills et al^ extended the spin ice label to include ma- 
terials based on corner-shared equilateral triangle mag- 
netic lattices, in showing that a local FM Ising model 
on the kagome lattice would have an even larger residual 
Pauling entropy: = iln(|). Ke et al^ experimen- 
tally demonstrated that the replacement of some mag- 
netic atoms from a pyrochlore lattice with nonmagnetic 
atoms indeed does produce a nonmonotonic variation of 
the residual entropy. One of us^ has recently shown that 
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TABLE I. Unit vectors ia denoting the local easy axis at each 
location within a unit cell for four inequivalent sites labelled 
by Q. The relative location of each site within the unit cell is 
given in terms of a crystal parameter u. 




FIG. 1. (Color online) The Trillium lattice. The direction of 
each spin is shown by an arrow (the direction of the spin at 
one of the four sites in the unit cell is pointing out of the page). 
Beyond the triangle a minimum of five bonds are required for 
the lattice to form a closed loop. 



in the disorder- free limit of such a replacement, hyper- 
kagome spin ice, one would expect the same enhanced 
Pauling entropy, = iln(|). 

In this paper, we report the first investigation^^ of spin 
ice on the trillium lattice. We show that a local FM Ising 
model with spins directed along the local < 111 > axis 
leads to a macroscopic ground state with a residual en- 
tropy [j^ « lii(f )] intermediate to those of kagome and 
pyrochlore spin ices, indicating that this is the newest 
member of the spin ice class. To motivate the experimen- 
tal search for trillium spin ice we calculate the magnetic 
susceptibility and neutron scattering structure factor for 
this model. We note the curious resemblance of the for- 
mer quantity to susceptibility measurements of EuPtSi^l, 
which remains disordered to very low temperatures and 
whose magnetic Eu^+ ions form a trillium lattice. We 
note that the neutron scattering structure factor is quite 
different from that of the AFM Heisenberg model^^ on 
the trillium lattice even in the cooperative paramagnetic 
phase and explain why this is expected to be a generic 
attribute of spin ice on corner-shared equilateral triangle 
lattices. 



where Sijocai = cTiaSa (no sum over a), a,a = ±1, (u) 
sums over nearest neighbors, J < is a FM coupling 
constant, and is a unit vector denoting the local easy 
axis of the trillium lattice as presented in Table H] As 
the dot product of any two of the nearest neighbor unit 
vectors ia is — |, this model is formally equivalent to an 
AFM Ising model on the trillium lattice. The trillium 
lattice is a simple cubic structure with a four site basis 
and a P2i3 symmetry featuring a 3D lattice of corner- 
shared equilateral triangles as pictured in Fig. [TJ 

It is interesting to note that none of the spin struc- 
tures corresponding to the ground state of the AFM Ising 
model are common to the ground state of the classical 
AFM Heisenberg model on any corner-shared equilateral 
triangle lattice in contrast to the case on the pyrochlore 
lattice. This means that the excellent agreemenlii^ be- 
tween large- studies of the AFM Heisenberg model and 
the local FM Ising model on the pyrochlore lattice is not 
expected to carry over to studies of corner-shared triangle 
lattices. 

In contrast to earlier studied spin ice lattices, the AFM 
Heisenberg model on the trillium lattice is known to or- 
der at low temperatures22ii2S despite the geometric frus- 
tration inherent to its structure. We will show below that 
this is not the case for trillium spin ice meaning that the 
FM local Ising model studied here is the first fully frus- 
trated model to be studied on this structure. The ground 
state of any corner-shared triangle structure can be sim- 
ply found by minimizing the Hamiltonian of Eq. [T] on a 
triangle and counting the energy per spin. For J < 0, 
this leads to a ground state configuration with two spins 
out(/in) and one in(/out) with a ground state energy of 
J per spin^. 



III. METHOD 



II. MODEL 

We consider a local FM Ising model on the trillium 
lattice with Hamiltonian, 

H = J ^S*;, local ■ Sjjocal, (1) 
{ij} 



Wc have carried out Monte Carlo simulations using the 
Metropolis algorithm with periodic boundary conditions 
on cubic lattices of side length L corresponding to AL^ 
spins. The bulk of our analysis corresponds to a choice 
of L = 6, which contains 864 spin sites, as our results 
appeared to reach the thermodynamic limit by L = 3. 
To ensure that we had reached the thermodynamic limit 
we considered systems as large as L = 18 in our calcu- 
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FIG. 2. (Color online) a) The acceptance rate of the spin 
flips as a function of the temperature. For L > 2, below 
about T = 0.2| J| the acceptance rate asymptotes to the L = 2 
exact result of 27.8% as discussed in the text, b) The only zero 
energy spin fluctuations allowed in the ground state feature 
spins whose neighbours all pair with opposite signs. Changing 
the direction of one of the sites changes three triangles from 
being mostly up to mostly down or vice versa as discussed in 
Sec. IV. B. 2. 



lation of the residual entropy of this spin system. At 
each temperature, 10000 Monte Carlo steps were used to 
equilibrate the system and a further 1000 were used to 
calculate the averages of physical quantities, where one 
Monte Carlo step was taken to on average attempt one 
update per site. Error bars as reported correspond to the 
standard deviation of our averages over four independent 
trials. Following the simulated annealing prescription of 
Kirkpatrick et ali^, for most of our analysis the temper- 
ature of the system was started at T = 20 J and reduced 
by 1% at each step. Exact results for the residual entropy 
were found by counting unique members of the ground 
state for the tractable periodic lattice sizes L — 1 and 
L = 2. 



A. Acceptance rate 

As shown in Fig. HJa) in contrast to studic&iS, of spin 
ice physics on the tetrahedron-based pyrochlorc lattice, 
at low temperatures on our triangle-based trillium lat- 
tice, the acceptance rate of the Monte Carlo simulations 
does not approach zero except in the case of L = 1 . This 
difference is easily understood in terms of the acceptable 
ground state spin configurations as shown in Fig. ^h). 
Such flippable sites do not exist on tetrahedron based 
spin ice lattices. Counting the total number of ground 
states for such a 3 triangle configuration which have a 
flippable spin gives 16 of 54 members of the ground state. 
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FIG. 3. (Color online) The evolution of the location of a 
particular flippable site in time. The path of spin flips which 
originated at site after a) 30, b) 109, and c) 191 updates. 
The numbers in the pictures are the number of updates it 
took for the path to reach that spin site. The number of 
updates includes updates where the flippable location does 
not change, instead the same site flips back and forth between 
its two possible directions. 



slightly above our low temperature asymptote. Such flip- 
pable configurations cannot exist for L = I. A more 
careful exact treatment for the 32 sites of i = 2 gives 
a total of 314874 ground state spin configurations. Of 
these 10075968 spin sites, we have numerically found that 
2797714 sites are flippable, yielding an exact result of 
27.8% of all sites arc flippable, for comparison with the 
i = 2 Monte Carlo asymptote of 27.9% ± 0.3%. As one 
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can see from Fig. ^a), there is no discernible difference 
in the acceptance rate for L = 2 and L = 6, indicating 
that we have quickly reached the thermodynamic limit. 

Lest one be tempted to suggest that for a given simula- 
tion one has roughly 73% of the spin sites exhibiting long 
range magnetic order amidst a background of static flip- 
pable spin sites, it is interesting to note that the location 
of the flippable sites is not fixed. Indeed, even at the low- 
est temperatures when the system can no longer access 
spin configurations outside the ground state spin mani- 
fold, after waiting a sufficiently long time every site on 
the lattice will flip. In this sense our model remains disor- 
dered (unlike pyrochlore spin ice which becomes stuck in 
a glassy stated) and able to equilibrate even to zero tem- 
perature. The dynamic process through which flippable 
sites move around the lattice can be illustrated by follow- 
ing the evolution of the location of one of the flippable 
sites as shown in Fig. [3] Here we see that it is possible 
to have neighbouring flippable sites. When a spin flip 
occurs at site one, there is some probability that a neigh- 
bouring spin at site two is now also a flippable site. When 
this is the case, the next flip could occur either at site 
one or site two. If the latter, the spin at site 1 may well 
no longer be flippable, and the location of the flippable 
site begins to migrate around the lattice. As seen in Fig. 
[3] this migration covers a reasonably large distance in a 
small number of updates in a manner reminiscent of a 
random walk or traveling salesman problem. 



IV. RESULTS 
A. Heat capacity 

Within a Monte Carlo simulation one computes the 
heat capacity as: 



a) 0.0002 



CviT) 



<E^> - <E>' 
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(2) 



where E is the energy of a particular spin conflguration 
on the lattice and T is the temperature of the system. 

Magnetically frustrated materials are distinguished by 
their lack of magnetic order to very low temperature de- 
spite the presence of strong spin-spin correlations. In 
particular, spin ice materials are known to exhibit a soft 
peak in the heat capacity indicating the onset of the ice 
rule at low temperatures. At high temperatures there 
is little energetic cost for flipping a spin in a disordered 
system. As more triangles begin to satisfy the ice rule, 
defects become more isolated forcing the system to pay 
an energetic cost to flip certain spins. Below the peak 
essentially all triangles are members of the ground state, 
with all accepted spin flip processes leaving the system 
in its ground state, thus paying no energetic cost. 

In Fig. m^a) we present the heat capacity per spin 
divided by temperature as a function of temperature for 
L — 6. Wc observe a characteristic soft peak in the 
heat capacity with the ice rule obeyed by the majority of 
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FIG. 4. (Color online) a) Heat capacity divided by temper- 
ature per spin versus temperature in units of the coupling 
constant J of the system at different temperatures for L = 6, 
b) Residual entropy per spin versus system size, L, for a cubic 
lattice with iL^ sites {L is the number of unit cells in each 
direction on the lattice). For L = 1 and 2 our Monte Carlo 
results agree with exact results. For L > 3 we appear to have 
reached the thermodynamic limit with an entropy 4% lower 
than than Pauling's approximation. Note that the region cor- 
responding to the Pauling entropy has been blown up to make 
our error bars visible, and the L — 1 result lies considerably 
higher. 



triangles below about -j, and completely obeyed below 
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B. Entropy 

At very high temperature one expects there to be es- 
sentially no correlations between the local Ising spins, 
or 2^ equally weighted spin conflgurations for an en- 
tropy per spin of ln(2). In a majority of materials, as 
the temperature approaches zero, in accordance with the 
third law of thermodynamics, the entropy of the system 
approaches zero. Spin ice materials are a notable ex- 
ception to this rule, as to very low temperatures they 
retain a large flnite residual entropy. One can quantify 
the amount of residual entropy at temperature T by sub- 
tracting the integrated weight under the ^ curve from 
the high temperature entropy of the system. 
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(3) 



Here and below we will refer to the entropy per spin for 
our model. In the T — > limit, we have a finite residual 
entropy and correspondingly macroscopic ground state if 
f(0) = ln(2)-/„°^~ 



1. Pauling estimate 

As noted above, the simplest approximation one can 
make to the residual entropy of a corner-shared trian- 
gle/tetrahedron lattice is to assume, as Pauling did in 
the context of water ice^S, that spins beyond the nearest 
neighboring triangle are uncorrelated. Then each spin 
site has two degrees of freedom constrained by only a 
fraction of the possible spin configurations on each trian- 
gle/tetrahedron being members of the ground state. For 
a triangle-based lattice this number is |. If the trian- 
gle assignments are uncorrelated then for each triangle 
added one adds a factor of | to the computation of the 
total states available to the system. On the hyperkagome 
and kagome lattices each spin belongs to two triangles, 
such that there are Nt = ^ triangles in total. On the 
trillium lattice each spin belongs to 3 triangles, so there 
are Nt = N triangles in total. The Pauling degeneracy 
of the ground state is then 2^(|)^' = 2^(|)^ = (|)^, 
hence the Pauling estimate for the entropy per spin is 
f = iln(|)^=ln(|). 



2. Exact 

While one might expect the Pauling estimate to pro- 
vide a reasonable approximation to the ground state de- 
generacy in the thermodynamic limit for uncorrelated tri- 
angles, it clearly has some limitations as the number of 
unit cells becomes small. As an example, if we consider 
the case L = 1, the Pauling estimate includes 4 triangles 
and 4 spins, leading to a non-integral number of distinct 
ground states, For this case it is not hard to show 
that all members of the true ground state have =4-1 



for two sites in the unit cell and ai, 



-1 for the other 



two. This describes (2) = 6 states, so the exact entropy 
per spin for the case L = 1 is = ;j hi(6), substantially 
higher than the Pauling estimate. 

The case L = 2 is still numerically tractable. If we 
designate the types of triangles available to the ground 
state by the number of "positive" (featuring aia = +1 
on two vertices, aia = — 1 on the other) and "negative" 
triangles (the converse), we see from Table lU that the 
geometry of the lattice only allows changes in the number 
of triangles of one type to vary by three triangles if they 
are to remain in the ground state. Counting all members, 
we find 314874 distinct states, much reduced from the 
Pauling estimate of approximately 431439.88 states. 



Positive Triangles 


Negative Triangles 


Number 


4 


28 


72 


7 


25 


1824 


10 


22 


19680 


13 


19 


76512 


16 


16 


118698 



TABLE II. The number of members of the L — 2 ground 
state featuring a given number of "positive" {aia = -f 1 on two 
vertices and aia = — 1 on the other) and "negative" {aia = — 1 
on two vertices, aia = -Hi on the other) triangles. The table 
is symmetric in positive o negative if continued. 



Although even the case L = 3 is no longer numeri- 
cally tractable, the necessity of changing the number of 
triangles of any type (positive or negative) by three con- 
tinues from one set of allowed triangle configurations to 
the next. This may indicate the presence of correlations 
not present in the Pauling estimate which could account 
for our results consistently falling somewhat below this 
approximation. The origin of this change of three trian- 
gles appears to be that a ground state configuration, in 
order to remain in the ground state with a single change, 
can only change a spin whose neighbors are all paired in 
opposite signs on their corresponding triangles as illus- 
trated in Fig. [Kb). 



3. Monte Carlo 

As shown in Fig. IHb) the entropy calculated by 
Monte Carlo simulations^ agrees well with the exact 
results available at small L, giving a residual entropy 
per spin at L = 1 of 0.44803 ± 0.00015 (c.f. with 
iln6 ~ 0.44794) and at L = 2 of 0.39564 ± 0.00010 
(c.f. with -Lln(314874) ~ 0.395623). As expected from 
Pauling's estimate, we find a large finite residual entropy 
as we approach the thermodynamic limit. For L ~ 3 and 
higher the entropy of trillium spin ice is consistent with 
0.3920 ± 0.0002, approximately 96% of the Pauling en- 
tropy (ln(|) ~ 0.40547), as one can see from the dashed 
horizontal line in Fig. |4l(b)^. 



C. Spin-spin corrlelations 

In order to calculate the static spin-spin correlations, 
the static structure factor (S) of the trillium lattice was 
found, 



id 

The dot product between all pairs of spins on the lat- 
tice was averaged in the same way that the heat capacity 
was, with 10000 Monte Carlo steps used to equilibrate at 
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each temperature, 1000 Monte Carlo steps used to find 
the average, and a temperature difference between each 
temperature step of 0.1%. The dot product between each 
pair of spins at a specific temperature was stored prior 
to multiphcation by the phase factor e*^ '-'''"''^-' where fi 
denotes the position of the spin within the lattice and 
hq the momentum transfer vector. The static magnetic 
susceptibility was found as 



i,3 



1. Magnetic susceptibility 

The temperature dependence of the inverse magnetic 
susceptibility {x~^) for L = 6 is presented in Fig. [5l[a). 
At high temperatures one sees (from the inset) that 
is to a very good approximation^^ a straight line with 
a best fit line = (0.991 ± 0.001)[-ji^]r - (0.49 ± 
0.01). Below about T = 3|Jcr^| one sees an upturn of 
relative to this line, with a weak shoulder around 
T = 'AJ£_\. whereupon x~^ approaches 0. 

2. EuPtSi 

In Fig. [Sfb) we show the inverse magnetic suscepti- 
bility of EuPtSi as measured by Adroja et al^. In this 
rare earth material, Eu^"*" is believed to lie on the tril- 
lium lattice (symmetry group of LalrSi, P2i3, and have 
nearest neighboring Ft and Si atoms lying along the lo- 
cal < 111 > axis. This suggests that as in other spin 
ice materials, single ion anisotropy may act to favour 
moments along the local < 111 > axis. As discussed 
below, it is not immediately clear why a Eu^+ (spin 
|, L = 0) state would experience strong crystal field 
effects pinning it to a local direction, although some 
evidence for an analogous favouring of a local trans- 
verse magnetization has been presented for equivalent 
Gd3+ ions in GdaTiaOy and GdaSnjOT^. Noting the 
qualitative similarity between Figs. El^a) and EJb), it 
is tempting to speculate that spin ice physics might 
be responsible for the lack of magnetic order seen to 
very low temperature in this material. To this end, we 
have fit the high temperature experimental data with 
X-\ = (0.1412 ± 0.0005) [^]r - (0.61 ± 0.09) [|^^. 

To convert this to SI units we need the unit cell length 
of EuPtSi: 6.436l2i. Noting that xsi = "^T^Xcgs 
Xsi = j^X7gsj ^-iid that there are 4 EuPtSi units per 
unit cell, we need to multiply our expression for x^^ in 

mol -u / 1 \ VNa _ 1 (6.536 xlQ-'^cm)^ (6.022 xlO^^mor^) ^ 
omu ^i-rr ' 4 47r 4 

3.194^j to arrive at dimensionless x^^ ^or comparison 
with our theory. Here V is the unit cell volume in cm"^ 
and Na Avogadro's number. This means Xsi EuPtSi — 
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FIG. 5. (Color online) a) Inverse magnetic susceptibility as 
a function of temperature in units of the coupling constant J 
times the square of the spin (cr^) from Monte Carlo simulation, 
b) Inverse molar susceptibility of EuPtSi from D.T. Adroja 
et al^. The inset of both graphs shows a linear fit of the 
high temperature points which was used to define the line 
that passes through the Curie temperature. 

(0.451 ± 2 X 10-3) [-^]T - (2.0 ± 0.3). From this informa- 
tion, it should be possible to extract an estimate of the 
coupling constant J, given a few assumptions. If Eu^+ 
takes a high spin configuration, one expects the spin to 
be s = ^, and orbital angular momentum to be zero, 
hence one would expect a Lande g factor of ^ = 2. If we 
ignore the discrepancy^ between the Monte Carlo and 
experimental x~^ intercepts, we can estimate based on 
the high temperature slopes of these lines an effective 
value of Jcr^ 0.455 ± 0.002^, where cr = |. 



3. Structure factor 

In Fig. Ela) and (b) we present the temperature evo- 
lution of the static structure factor in the hhk and hkO 
planes respectively. One sees that correlations gradually 
increase in strength until entering the spin ice regime. 
For comparison with previous work we have set the pa- 
rameter u = 0.138, the value relevant for the Mn sites 
of MnSi, as the value of u for EuPtSi is not known. 
Qualitatively correlations are seen to be quite different 
in distribution from those seen for the AFM Hcisenberg 




FIG. 6. Relative intensity of the structure factor in the a) hhk and b) hkO plane. Temperatures decrease from top left (10.266J), 
to top right (1.787 J), to bottom left (0.709J), to bottom right (0.0642J). The axes range from to ^ where a is the 
lattice constant with k on the vertical axis and h on the horizontal. The structure factor was calculated using it = 0.138 as 
discussed in the text. The relative strength of the correlations is shown on the right of the graphs with white indicating areas 
of strong correlation and black indicating areas of weak correlation. Although the color scale remains the same for each image, 
the number of divisions is 7 for the lowest temperature, 10 for 0.709J, 50 for 1.787J, and 90 for 10.266J. 



model on the trillium lattice as shown in Fig. 3 of Isakov 
et al^. In particular, one sees a lot of weight inside the 
box defined by { — ^, ^} in each coordinate in the spin ice 
state, while there is very little if any weight for the AFM 
model. At first the difference between the spin-spin cor- 
relations of these models may appear surprising as one is 
used to the very good agreement of spin correlations seen 
on the pyrochlore lattice for a local FM Ising model with 
mean field calculations of an AFM Hcisenbcrg modeU^. 
However, as explained above this is a simple difference 
between corner-shared triangle lattices and corner-shared 
tetrahedral lattices - that elements of the ground state 
manifold of an AFM Ising model are not common to the 
ground state of an AFM Heisenberg model on corner- 
shared triangles. Put simply, the sum of the spins on 
a triangle in an AFM Ising model cannot be 0, which is 
the condition for the ground state of the AFM Heisenberg 
model on a triangle. 



V. DISCUSSION 

In light of the recent interest in the creation of de- 
confined magnetic monopoles on spin ice lattices, it is 
interesting to consider whether or not deconfined mag- 
netic charges could be present on any of the three corner- 
shared equilateral triangle spin ice lattices: kagome, hy- 
perkagome, and trillium. Certainly, on each of these 
lattices, the nature of the ground state dictates that 
the coarse grained V • B = condition responsible 
for the "bowties" seen in polarized neutron scattering 



experiments^ of Ho2Ti207 cannot hold true for all mem- 
bers of a ground state as in its most basic version it re- 
lies on the sum of the spins on each triangle/tetrahedron 
vanishing. While this can in principle occur for AFM 
Heisenberg spins on a triangle based lattico22,, the ground 
state of the FM local Ising model cannot satisfy this con- 
dition. One might imagine a more coarse-grained sum 
over several spins could be satisfied for some members 
of an Ising ground state, with total spin 0, but it seems 
highly unlikely that such a coarse graining could capture 
all members. As such, one might not expect to see the 
"bowtie" structures common to the neutron scattering 
structure factor of AFM Heisenberg models on the py- 
rochlore, kagome, and hyperkagome lattices in the spin 
ice analogs on kagome and hyperkagome lattices. Given 
this, one might find it surprising that dipolar spin ice ap- 
parently does exist on both of these lattices^i^. It is an 
interesting open question whether dipolar spin ice exists 
on the trillium lattice, and whether or not the ground 
state excitations on any of these corner-shared triangle 
lattices can be thought of in the language of deconfined 
magnetic charges. 

On the experimental side, it is worth noting that the 
trillium lattice is a reasonably common magnetic sub- 
lattice featuring an itinerant FM (MnSi), a Kondo in- 
sulator (FeSi), and several large moment Eu compounds 
(EuPtSi, EuPdSi, EuPtGe, and EuIrP) not to mention 
a polar molecule CO which exhibits a residual entropy 
slightly smaller than In 2 per spiniS. Of the Eu com- 
pounds magnetic order has not been seen to below 4.2K 
in the first two, to below 2K in the third, with only the 
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fourth seen to ferromagnetically order—. To date the first 
three materials have been classified as paramagnets, but 
it is our belief the possibility that they are cooperative 
paramagnets has not been considered. In this context, 
it may not have been realized that the magnetic sublat- 
tice of these materials is actually a frustrating lattice. 
The presence of a cooperative paramagnetic phase would 
help explain why each of these materials appears to have 
a small yet positive Curie temperature. 

While we are struck by the qualitative (and perhaps 
quantitative) agreement between our model and mag- 
netic susceptibility data of EuPtSi, it is hard to under- 
stand why a local FM Ising model would result from a 
detailed investigation of this material, given that a half- 
filled spin shell should be spherically symmetric if the 
high spin state is adopted. We mentioned that the Eu 
atom lies between Pt and Si atoms along the local < 111 > 
symmetry axis of the crystal, however there are also 2 
nearby triangles of Pt and Si atoms centred along this 
axis which bring the coordination number of the Eu sites 
up to 20 when the nearest neighbor Eu atoms are consid- 
ered. Further, the related material EuPtGe is metallic to 
quite low temperature^ (the resistivity of EuPtSi has not 
yet been measured), suggesting that any magnetic cor- 
relations that do arise in EuPtSi (if it is metallic) might 
have a Ruderman-Kittel-Kasuya-Yoshida nature^. Pre- 
vious work by one of us^ has shown that the extended 
Heisenberg model with Ji FM and J2, J3 AFM has re- 
gions of momentum space which remain nearly degener- 
ate to low temperatures, suggesting that it may be pos- 
sible for magnetic frustration to play a role in the sup- 
pression of the ordering temperature of EuPtSi even if it 
is not a true realization of the spin-ice physics we have 
described in this work. 

Lastly, we should mention the possibility that artifi- 
cial versions of spin ice on the trillium lattice could be 
created if it turns out that none of the naturally occur- 
ring trillium ice candidates provide an immediate physi- 
cal realization of this model. When kagome spin ice was 
proposed. Wills et aZ.— could hardly have predicted that 
the first realizatioi*^ of their model would take the form 
of lithographically etched mesoscopic islands, whose only 
magnetic interactions were dipolar in nature. While it is 
hard to imagine lithographically etching a 3D lattice it 



is not inconceivable that such technologies could evolve. 
Likewise, advances in optical lattices have experimental- 
ists claiming^ that in the near future model spin systems 
may be able to be built to order, perhaps yielding an al- 
ternate path to the realization of the interesting physics 
of trillium spin ice. 



VI. CONCLUSIONS 

In this paper we have predicted the existence of spin 
ice on the trillium lattice. We have shown that a lo- 
cal FM Ising model on the trillium lattice has a macro- 
scopic residual entropy of 0.3920 ± 0.0002 per spin. We 
have calculated the heat capacity of this model and have 
seen that it takes the characteristic shape of spin ice on 
other lattices, with no signs of magnetic order to the low- 
est temperatures. The magnetic susceptibility at high 
temperatures appears to follow a Curie Weiss law with 
a positive Curie temperature before smoothly deviating 
toward = 0, exhibiting an intriguing weak shoulder 
at around ||J(T^|. We have identified a material whose 
magnetic susceptibility appears to show some of these 
features. The neutron scattering structure factor of this 
model has been seen to noticeably differ from the AFM 
Heisenberg model on this lattice, a feature we expect to 
be generic of corner-shared triangle lattices. It remains 
an interesting open problem whether or not the dynam- 
ics of spin excitations on this and other triangle-based 
spin ice materials can be well-described in terms of free 
magnetic charges, and whether dipolar spin ice exists on 
the trillium lattice as it does for other spin ice lattices 
investigated to date. 
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